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TECHNICAL NOTE 5296 

SEPARATION, sWeILITY, AND OTHER PROPERTIES OF 
COMPRESSIBLE LAMINAR BOUNDARY LAYER WUM 
■ PRESSURE GRADIENT AND HEAT TRANSFER 
By Morris Morduchow and Richard G. Grape 


SUMMARY 


A theoretical study is made of the effect of pressure gradient, 
wall ten5>erature, and tfech number on 1.am1 nar boundary- layer chEiracter- 
Istics and, in paarfclcular, on the sMn-frlction and heat-transfer 
coefficients, on the separation point in an adverse pressure gradient, 
on the ira.ll temperature required for conqplete stabilization of the 
laminar boundary layer, and on the mlnlTnum critical Reynolds number for 
laminar stability. The Prandtl number is assumed to be unity and the 
coefficient of viscosity is assumed to be proportional to the tempera- 
ture, irLth a factor arising from the Sutherland relation. A simple and 
accurate method of locating the separation point in a compressible flow 
with heat transfer is developed. Numerical examples to illustrate the 
results in detail are given throu^out. ' 


INTRODUCTION 


The piirpose of the present investigation is to determine theoreti- 
cally the nature of the laminar boundary layer in ccmpressible flow with 
heat. transfer and pressure gradient. paarbicular, the effect of 
pressure gradient (favorable and adverse), wall tenperature, and Mach 
number on the boundary-layer characteristics are investigated. Such an 
investigation has already been made in reference 1 on the basis of 
fourth-degree velocity and stagnation-enthalpy profiles, in conjunction 
with two different boundary-layer thicknesses (a dynamical and a thermal 
boundaiy-layer thickness) . In contrast with the present study, the 
effect of normal fl\ild Injection at the wall was included in reference 1 , 
but the stability of the boundary layer was not investigated therein. 

The present investigation is based on the more accurate method of calcu- 
lating boundary-layer properties as develcQped in reference 2, where 
sixth-degree velocity and seventh-degree stagnation-enthalpy profiles 
are applied in conjunction with a single boundary-layer thickness, with 
the thermal boimdary- layer thickness replaced by an additional param- 
eter b2_ in the stagnation-enthalpy profile. 
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Iq the first section of this study, skin-friction and heat-transfer 
characteristics are investigated, especially with respect to the effect 
of wall temperature, pressure gradient, and, Mach number. Ihe analysis 
is carried out first in general terms and then Illustrated hy a numeri- 
cal example for the supersonic flow over a thin airfoil, for which the 
pressxire gradient is everywhere favorable (negative) . In the next 
section, a simple and ordinarily sufficiently accurate method of calcu- 
lating the separation point in an adverse pressure gradient in subsonic 
or supersonic flow^ over a wall at a given uniform tenperature is devel- 
oped. The method, based on the special use of a seventh-degree profile 
to satisfy an additional boundary condition at the separation point 
(first suggested in ref. 5 ) j is essentially an extension of the method 
developed for zero heat transfer in reference It, where a numerical 
exanple indicated excellent agreement with exact results. With this 
method, the effect of wall temperature and Mach number on the separation 
point is investigated- The results are then illustrated by a numerical 
exanple for flow with a linearly decreasing velocily outside of the 
boundary layer. A second exanple treats the conditions (involving the 
wall tenperature) under which laminar separation will take place imme- 
diately behind a stagnation flow abnq)tly followed by an adverse 
pressure gradient. 

In the final section, the stability of the laminar boundary layer 
over a thin biconvex airfoil in stapersonic flow is investigated. For 
this purpose, the wall tenperature required for complete stabilization 
of the flow (infinite minimum critical Reynolds number) is calculated 
for several Ifech numbers at different stations along the airfoil. In 
addition, the minimum critical Reynolds number ' is determined at a given 
station for various wall temperatures and Mach numbers. ConpariBon is 
made with results for flow over a flat plate in order to demonstrate 
the effect of a pressxire gradient, in addition to that of wall tenper- 
atxire, on the stability characteristics. The stability calcxilations 
here are based on the well-known two-dimensional criteria developed by 
Lses and Lin (refs. 5 and 6 ) for compressible flow. Since the validity 
of these criteria appears at present to be in doxjbt for hi^ Mach 
numbers (e-g., ref. 'j) , the present calcxilations have been restricted 
to Mach nxmibers not above 5 * The present resxilts on stability charac- 
teristics in conpresslble flow with heat transfer and pressxire greidlent 
may, in a sense, be considered as an extension of the resxilts already 
obtained for flow over a flat plate (zero pressxire gradient) with heat 
transfer (refs. 8 to ll) and for flow with a presBxire gradient but with- 
out heat transfer at the wall (refs. 12 and h). 

The analysis in the present investigation is based on the assxmption 
of a xmiform wall tenperature, a Prandtl nxmiber of xmity, and a linear 


-*^ock-wave interaction with the boxmdary layer, however, is not 
considered. 
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temperature-viscosity relation (as in ref. l) - As in reference 2, how- 
ever, a factor C (first suggested and applied in ref. 15 ) has heen 
introduced in the temperature-viscosity relation to account, at least 
at the wall, for the Sutherland viscosity law. This factor, as will he 
seen, has an influence on the skin-friction and heat-transfer coeffi- 
cients, aJLthou^ not on the separation point. One sin^jle, hut very 
approximate, means of correcting the Nusselt numhers to he obtained 
here in cases of a Prandtl number different from unity would he to 
multiply the Nusselt numbers for unit Prandtl number hy a power of the 
actual Prandtl number, this power being roughly equal to I /5 (cf., e.g. 
refs. 14- and I 5 ). 

This investigation was conducted at the Polytechnic Institute of 
Brooklyn Aeronautical Laboratories under the sponsorship and with the 
financial eisslstance of the Ifa,tio 3 jfe.l Advisory Committee for Aeronautics. 
The authors hereby e 3 q>ress their thanks to Professors P. A. LLhhy and 
M. Bloom for their helpful discussions and to Ifr. Richard P. Shaw for 
his aid in the calculations. 
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ratio of stagnation enthalpy at wall to stagnation 
enthalpy at outer edge of boundaiy layer, 

for a Prandtl number of 1 , is ratio of actual 

w a l l temperature to equilibrium wall temper- 

ature Tg for zero heat transfer (see also eq.. (l 5 )) 

stagnation enthalpy, defined as quantity (u^/2) + c^T 
constants (see eqs. (28) and (29)) 
thermal, conductivity 

coefficients in velociiy distribution over thin airfoil 
(eqs. (it)) 

characteristic lengthy chord length for airfoil of 
figure 1 

Mach number v 

Ifuaselt number 
Reynolds number, vi^L/Voo 
Reynolds nimiber, vtjjLyiq) 

minimum critical Reynolds mmiber, "based on conditions 
at point b immediately behind shock wave at leading 
edge of airfoil in figure 1 , (ubL/Vb)cr 

minimum critical Reynolds number, based on remote free- 
stream conditions in siq)ersonic flow over thin airfoil 
of figure 1 , (UooL/v„) 

proportional to ratio of local skin friction to local 
Nusselt number (eq. (16)) 

constant in Sutherland vlscosity-teii5»erature relation, 
216 ° R for air 

absolute tengperature 

equilibrium wall temperature for zero heat tremsfer j 
for a Prandtl number of unity, Tg = - 


t 


variable defined by equation (5) 
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velociijy in x-dlrectlon 

N 

coordinates parallel and. normal to surface, 
respectively 

ratio of specific heats, Cp^c^, lA for air 
density 

■boundary- layer thickness in xt-plane 
geometric slope of thin airfoil profile 


[i, coefficient of viscosity 

V kinematic viscosity 

I dimensionless distance along wall, x/L 

I’ = x‘/L (see fig. 1 ) 

T dimensionless variable, tyS-t 

q>i constant defined hy equation (8) 


Subscripts: 


s,sep 


constant used Instead of cp^ in determining separation 
point (see eq. (S^t)) 


values at point outside of boundary layer immediately 
behind shock wave at leading edge of airfoil 

critical 

favorable presstnre gradient 
values at wall 

values at or used for determining separation point 
(see eqs, (25) and (21)) 

local values at outer edge of boundary layer 
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00 values at suitable reference point outside of "boundary 

layerj for airfoil in supersonic flow, values in 
undisturbed (remote) free stream 

' differentiation with respect to g 


BASIC EQimnnioNs 


The "basic equations for the calculation of the laminar-houndary- 
layer characteristics in coii 5 )res 8 ihlje flow with an axial pressinre 
gradient over a wain, at a uniform wa l l ten^ierature Tq have been 
developed in reference 2 with the assumptions that the Prandtl nimrber 
is Tinity and that the coefficient of viscosity varies with tenperature 
according to the relation 


= c(T/o;„) 

(1) 

where 


C = (V^co)^/^(2^ + S)/(To + S) 

(2) 

These eq.uations will be repeated here for convenient reference. 
The variable t replaces the physical normal coordinate y according 
to the relation 

y = r (T/Ti)dt 

(5) 

The velocity and stagnation-enthalpy profiles are given, 
■by 

respectively. 

i = ^2 t - + 6x5 _ 2t^) + - loA + 

10x5 - 3x^ + 

+ 10x5 _ 20x^ + 15x5 - l4-x^^ 

(4) 

^ + (1 - Gi)(55t^ - 81^x5 + 70 x^ - 20 t'^ ) 

+ 

bi(x - 20 x^ + 45t5 - 56x^ + lOr"^) 

(5) 
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where 




(6) 


and the dimensionless boundary-layer thickness A(|) in the xt-plane 
can he calculated from the eq.uation 


ropj'""’ 


[(27-1)/ (r-i: 

1] - 

1 



(7+1)/ (7-1)] 

- (‘Pl/^ 

) 


( 7 ) 


Here and are constants defined as 

= 0.1093 + 0-002Ua2 - 0-000622 s2^ + 0.0004l2^bj^a2/Gi) - 
0. 0000095 (^ia2/li)^ - 0.000155 (ibiag^/Gi) 

= 0-5Gl + O.OCA-58 + O.O 232 S 2 - 0.00124s2^ + 

ti|o.0905 + (agy^OGi) (0.0858 - 0.00458ag) 

where S2 and bj^ denote constant average values of and b^Cl) 

over the entire flow. 


( 8 ) 


Equation (7) is a solution of the ordinary differential equation 

r - Fi)j I 


(fl/2)7v’ + x|Pi(pi'/Pi) + (uiyui)[^5>l + 


( 9 ) 


Tlip coefficient b3^(^) is a function of 8 ^{i) and is to be calculated 

as the solution of a quadratic equation (eq. (68) of ref. 2) after A(S) 
and thence a^i^) have been determined. It will usually be found that 

an approximate value for b^^ is 


bi « 2(1 - Gi) 


(10) 
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(Eq. (10) Is exactly valid for flow over a flat plate,) The temperature 
dlstrihution Ti/1^ and the Ifech ntmiber at the local outer edge of 

the houndary layer will he related to the velocity dlstrihution 
there according to the relations 

TlA„ = 1 + [(7 - l)/2]t^2[l - (11) 

= K/^oo)V(Ti/I;>)-'^ (12) 

The wall-te3i5)erature ratio '^ q/%o is related to the parameter = Tq/Tq 
hy the relation 

= Gi(i + (15) 

The local skin-friction and heat-transfer (Nusselt number) coeffi- 
cients follow, respectively, fran. the eqmtions 

Cf J = - (=2/5) - (‘> 1 “ 2 / 600 i)] (c/lfX)(T3/3i.)(ui/l„)V^/2 

(W 

SKIN-FRICTION AND HEAT-TRANSFER CHARACTERISTICS 


From the foregoing equations, several general conclusions on the 
effect of wall tenq)erature, pressuore gradient, and Lfech number on the 
skin-friction and Nusselt number can he derived. 


General Dnpll cations 

The effect of wall tengierature on the skin-friction and heat-transfer 
coefficients arises essentially from two different 80\urces: (a) Ihe 

conditions of dynamic equilibrium, as defined mathematically hy the basic 
differential equations and their solutions, and (h) the variation of the 
viscosity coefficient with tenqperature, as defined mathematically hy the 
constant C. 
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With respect to (a) , the effect of wall ten 5 )erature will depend on 
the natinre (positive or negative, i-e., adverse or favorable) of the 
pressure graxilent. This follows from the fact that in the governing 
momentum ordinary differential equation (eq. ( 9 )) the quantity Gi 
appears primarily in a form multiplied hy the velocity gradient Uj* 

in the flow outside of the hoimdary layer. This is also the case in 
the coefficient 02 (eq. (6)). By bringing the U]_' term in equa- 
tion ( 9 ) to the right side and observing the expression for according 

to equation (8), it can be seen, in fact, that for uj^' > 0 (favorable 

pressure gradient) (T^C) ', and hence {'h/C) , will tend to be increased 
by a decrease in the wall-tenperature parameter This result, in 

conjunction with equations (iJt-), (l5)^ ( 6 ), and (lO), indicates that, 
without the effect of C, lowering the wall temperature tends to diminish 
the local Nusselt number and, especially, the local skin friction in a 
favorable pressure gradient and to increase the Nusselt number and, 
especially, the skin friction in an adverse pressure gradient. This 
conclusion is in accord with that also derived in reference 1. It must 
be noted, however, that for a cooled wall, where G]^< 1 , the foregoing 
effect of wall temperature on 82 due to a pressure gradient will 
ordinarily be greater than that on A/C, since it will be found from 
the general solutions ( 7 ) and (8) that A/C is then not very sensitive 
to changes in G]^. 

From equations (l4), (15)^ and ( 7 ) it follows that both the skin 
friction and Nusselt numiber will be proportional to fU". The effect of 
wall tenperature arising from the viscosity coefficient and determined 
by the constant C is, contrary to the dynamical effect just discussed, 
independent of the pressure gra^ent. From equation (2) it follows that, 
if To/Ha is diminished, then C and hence Cf ^ and Nu decrease 

when To/%o< (or Tq < 2l6° R) , while C and hence Cf and Nu 

increase when Tq/ 5^ > (or Tq > 2l6° R) . The latter is expected 

to be the case in practice. 53ius, depending on the nature of the 
pressvire gradient and possibly on the magnitude of the wall temperature, 
the dynamical and viscosity effects of wall temperature may tend either 
to magnify or partly to cancel each other. More specifically, unless 
the wall is unusually cold, that is, unless Tq < S, these effects will 
tend to oppose each other in a favorable pressure gradient and to magnify 
each other in an eidverse pressure gradient. 

p - - , — . 

'^According to the equations developed here, this conclusion may not 
be quite valid in the immediate vicinity of a sharp leading edge where 
the pressure-gradient effect is_relatively uniuportant and the small 
effect of wall temperature on F]_ may actually predominate. (Cf. table II 
in conjunction ■td.th the niimerlcal example discussed subsequently.) In 
this vicinity, hcnfever, the dynamical effect of a uniform wall tempera- 
ture will be found to be negligible. 
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It is significant to note that, since the velocity gradient 
appears in the equations (eqs. ( 6 ) and (9)) always multiplied hy the 
wall-teii5)erature ratio G]_, it can he inferred that a lowering of the 
wall ten5)erature has a tendency to diminish the direct effect of a given 
pressure gradient, that is, the effect of u^’ as such, on the houndary- 
layer properties. A clear illustration of this will he seen subsequently 
in the analysis of 1 amlnar separation (cf . also ref. l) It must he 
observed, however, that the effect of a pressure gradient also appears 
indirectly , namely, in the variation of u^yu^o with |. 

For l&ch numbers above 1 , in fact, the Ti/Too terms in Tv (eq. (7)) 
may become particularly In^jortant, so that in such a case the net effect 
of the pressvtre gradient may actually he increased hy a lowering of the 
wall temperature. This will he clearly illustrated hy the suhseqixent 
numerical example. 

From equations ( 14 ) and (15) it follows that the ratio of local skin- 
friction coefficient to the local Russelt nuniber cein he expressed as 


r 



4(1 - Gi) 
^1 


1 



(16) 


Equation (16) is valid along the entire flow. For flow over a flat plate 
(zero pressure gradient: ~ ^ ~ O) , equation (16) implies 

r = 2 . For flow with a pressiure gradient, however, this simple relation 
is seen to he no longer valid. Since, ordinarily, 4(1 - Gi) /h^ » 2 , 
eqxiations (16) and (6) imply r > 2 along the flow in a favorable 
pressiire gradient (u^* > 0) and r < 2 in an adverse pressure gra- 
dient (u^' < 0 ) . Thus the ratio r tends to he increased hy a negative 
(favorable) pressure gradient and decreased hy a positive (adverse) 
pressure gradient . 4 From equations (6) and (16) it is seen that, for a 


5 Further illustrations of this conclusion can he found in numerical 
examples of reference I6, which are based on a small-perturhation method. 

4 Equatlon (l 4 ) and hence eqviation (16) may he inaccurate at, and 
hence in the immediate vlcini~ty of, the separation point in an adverse 
pressure gradient, since, according to the criterion of separation 
(eq. (20)) developed in the succeeding section "Location of Separation 
Point," Cf^ as given hy equation (l 4 ) will not vanish exactly at the 

separation point. This is due to the lise of sixth- instead of seventh- 
degree velocity profiles here. However, equations (l 4 ) and (16) should, 
for practi-cal purposes, he adequate to yield the distribution of skin- 
friction along the flow. When the primary interest is in the location 
of the sepaiation point, the method described in the section "Location 
of Separation Point" shotOd he used. 
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given pressure gradient and Mach number, the lower the wall ten^ierature 
the closer will r be to its value for flcnf without a pressure graidient. 
This Illustrates the diminution of the direct effect of a pressure gra- 
dient by cooling of the wall. 

Jl’can eqxxatlon ( 7) ^ as has already been noted, it will be found that 
in the presence of a pressure gradient 7 \/C may be appreciably affected 
by the Mach nximber because of the values of Ti/Tix>(l)- Consequently it 
can be inferred, in view of equations ( 1 ^ 4 -) and (15)^ that a pressure gra- 
dient will in general tend to enhance the effect of Mach number on both 
the heat-transfer and skin-friction coefficients. Since the values of 
Ti/T^(|) will depend on the distribution of the velocity 

along the flow (eq. (U)), this effect of Jfech number will, in fact, 
depend on the nature of the pressure gradient. For a favorable pressure 
gradient, for exangile, one for which (with proper choice of Uj„) 
ui/uoo > 1 and hence Ti/Q^ < 1 , an increase of Mach number can be 
expected to Increase 'k/C and hence, according to equations (l 4 ) and 
(15), to decrease both EFusselt number and the skin friction. The oppo- 
site effect will tend to occur in an adverse pressirre gradient. 

It should be further noted that for a given ratio G]_ of wall 

tenq)erature to equilibrium tenq)erature a lyb.ch number effect, independent 
of the pressvtre gradient, also appears in the viscosity-temperature 
factor C according to equations ( 2 ) and (15) . If Tq > S, then for a 
fixed Gi an increase of Mach number will diminish C and hence 

will tend, as far as C is concerned, to diminish both the skin-friction 
and heat-transfer coefficients in preportion to l/^. Thus, in a favor- 
able pressucre gradient, the dynamical (i.e., Tp/H^) and the viscosity 
(i.e., C) effects of Mach nunber will tend to amplify each other, while 
in an adverse pressutre gradient they will tend to oppose each other. 

It may be worth while to note here that, in view of the fact that 
A, and hence 6 ^, will ordinarily be only little affected by the wall 
temperature, equation (5) implies that cooling of the wall will in 
general tend to diminish the physical bo^Iadary-layer thickness. However, 
for a given value of Gp, the boundary-layer thickness will tend to 
Increase with Mach nunber. 


numerlceil Example 

In order to illustrate the foregoing general conclusions the 
boundary layer in the supersonic flow over a thin biconvex circular-arc 
airfoil of thickness ratio O.CA- (fig- l) at zero angle of attack was 
calculated by means of the equations given here. 5 For this case, the 

^ 1^1(1) # however, was calculated by means of equation ( 68 ) of refer- 
ence 2 and not by the approximate equation (lO) given here. 



12 


MCA TN 5296 


velocity distribution outside of the boundary layer can be exj)re 6 sed 
as (ref. 17) 

^ 1 /^ = 1 - kin - kgn^ 


1^1 = (^.2 - 1 )"^/^ 


(IT) 


k2 = 





where n is the slope of the airfoil at any point along its surface. 
In equations (I7), the subscript <» refers specifically to the remote 
free-stream conditions. 


It is in this case convenient, especially for purposes of compar- 
ison with flow over a flat plate, to use the point (to be designated by 
subscript b) Immediately behind the shock wave at the leading edge as 
the reference point. 6 For this purpose, it may be noted that 

iil/ub = 

where ui/u^o is given by equations (17)^ •while ujj/u^, is ob-fcained from 
eqmtlons (17) by e'valuating the slope n leading edge. The 

free-stream I^fech number can be e:iqxreased in terms of the Mach 

number at point b by measas of the relations 

Since in equations (l) to (16) -the subscript » refers to any sui'table 
reference point outside of -tiie boundary layer and since in siqpersonic 
flow it must refer, in "these equations, to a point behind a leading-edge 
shock wave, it -will, for this exangple, be taken to deno-te the point b. 

It -will be assumed in -this exanple -that S/Tj, = 0.14-16. T It should be 
noted -that the pressure gradient here is entirely favorable (u^* > 0 ) . 

6 in this manner, when con5)arison is made wi-th flow over a flat pla-te, 
the effect of 'the pressiire gradient over 'the airfoil -will be re'tained, 
but "the effect of •the leading-edge shock 'wave ■will be essentially e li minated. 

Tin "this example, -therefore, "the "tempera-ture outside of the 

boundary layer immedla-tely behind -the lesuiing-edge shock wave, ins"tead of 
•the remo-te free-s-tream -tempera-ture 5 ^, is considered fixed, while "the 
Mach number behind -the shock -varies. 


(18) 
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Figures 2 and 5 show the distrihutlon of local Nusselt ntmiber and 
skin friction for this example, with = I .5 ancL 5-0 and with 
Gl = 0.3^ 0.5^ and 0.7- results for a flat plate are also shown 

for comparison. (For = 1.5^ the value of 02 was taken to he 
02 = -I. 62 OG 1 , while, for = 3-Oj the value ap = -I. 55 OG 1 was used.) 
Tables I, H, lU, and IV give the calciolated values of C, 7\/C, 

respectively, for the various values of and 

G]^. In accordance with the general conclusions developed here, it should 

he expected that, for this exan^ile of a favorable pressure gradient, the 
dynamic effect of wall temperature (i.e., the effect of wall tender ature 
if C were fixed) shoiald he such that lowering the wall temperature 
diminishes both the ITusselt number and, especially, the skin friction. 
However, since with the values assumed in this example T 0 /T 5 > S/T 5 in 
all the cases, the viscosity (or C) effect of wall temperature will here 
be opposed to the dynamical effect (cf . also tables I and U) . Figure 2 
indicates that the Nusselt number here increases as the wall temperature 
is lowered for a given Mach number and hence that the effect of wall 
ten^jerature arising from the viscosity- temperature relation is the pre- 
dominant effect here. For the skin friction, however, figure 3 indicates 
that except near the 3jeading edge (5’ < 0.2) the dynamical effect of the 
wall ten 5 )erature is predominant for = 1 . 5 ^ since here Cf^ diminishes 

as the wall is cooled. For = 3*0, on the other hand, where Tq/Tjj 
is greater for a given value of (cf. eq,. (13))^ so iliat the C-effect 

becomes more important than for = 1 . 5 ^ it is seen that the dynamical 
or pressure-gradient effect of wall temperature predominates only slightly 
and only after a considerable distance (^’ > 0 . 6 ) downstream of the 
leading edge. 

Figure 3 indicates that, for = 1.5^ lowering the wall ten 5 >era- 
ture brings the skin-friction curves for the airfoil closer to those for 
a flat plate. This Illustrates the general conclusion previously reached 
regarding the diminution of the direct effect of a pressure gradient by 
a lowering of the wa3JL temperature. However, for = 3-0, figure 3 
indicates that the skin-friction curves for the airfoil will now be 
further removed from those for a flat plate when the wall temperature is 
lowered. This is essentially due, as previously intimated, to the 
Increased effect, at this higher Mach number, of the temperature, distri- 
bution T 2 yTb(^)(<l) outside of the boundary layer along the airfoil, 
which tends to increase ?y/C and hence diminish Cf^ (as well as Nu) 

in comparison with the value over a flat plate (where T^/Tb = l) • The 
direct effect of the pressure gradient given by the 02 term in equa- 
tion (l 4 ), however, is to increase the skin-friction coefficient over 
that for a flat plate, but at the higher Mach number this effect is not 
so great as the indirect effect of the pressure gradient due to T]^/Ttj(|). 
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The decarease of local Niisselt number for "both = 1.5 and 
= 5*0 (fig. 2 ) and the decrease of the local skin-friction coeffi- 
cient for = 5-0 (fig. 3 ) along the airfoil downstream of the 
leading edge are due primarily to the decrease, in this region of nega- 
tive pressure gradient, of the local ten 5 )erature outside of the 

boundary layer (which increases l/C according to eq.. ( 7 )). This, in 
fact, illustrates the enhanced effect of Mach nuuiber on the heat- 
transfer and skin-friction coefficients due to the pressure gradient. 

Table HI indicates that the values of b]^ as calculated from 
equations (6ja) and (68) of reference 2 remain fairly close to the 
value 2(1 - Gp) (cf. eq. (lO)), although in the present favorable 
pressure gradient (a 2 < O) they are everywhere less than or equal (at 
the leading edge) to this value. 


LOCATION OF SEPARATION POINT 


In a region of adverse pressxure gradient (negative u^') there is 
a possibility of laminar separation, which occurs where (hu/by) ^ = 0. 

A fairly accurate and simple method’ of calculating the separation point 
in compressible flows with zero heat transfer was developed in refer- 
ence 1+ and was baned in part on the use of an additional boundary con- 
dition (first suggested in ref. 5 ) at "tiie wall necessarily satisfied at 
the separation point by an exact solution of the partial differential 
equations. For zero heat transfer, the separation point as a function 
of Mach number calculated by this method was found to agree very well 
with nuiffisrical solutions in reference 18. (See ref. 4 for details; 
also, see table V.) In the present section this method will be general- 
ized for any given uniform wall -temperature ratio Gp = To/Te* (For 
zero heat transfer, as in ref. 4, Gp = 1.) 


Ifethod of Calculating Separation Point 

Py differentiating the momentxmi partial differential equation of 
the laminar boimdary layer, it can be shown (see the appendix), under 
the present assimptlons of a Prandtl number of 1 and a linear viscosity- 
temperature relation, that at the separation point 

= 0 ( 19 ) 

■^fith, as well as id-thout, heat transfer at the wall. 
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A seventh-degree° velociiy profile in t can now be chosen to 
satisfy condition (I9) in addition to the boundary conditions satisfied 
by the sixth-degree profiles (eq- (4)) on which the preceding analysis 
has been based. From this profile (cf. appendix) it follows that separa- 
tion will occur where haa the value (denoted by aps) 


^2s 


Gi + (2/l5)bi 


( 20 ) 


Moreover, is still given by equation (6). From equations (6) 

and (20) it follows that the value Aggp of A at the separation point 
will in general be 


^ep “ 


(^1/^) 


“(2-7)/ (7-1) 


K'/^)(l + Gl + (2/15)151 


( 21 ) 


By inserting the approximate equation (lo) for bj^ into equation (21),^ 
it is found that 


Isep = -105c 


(Tl/To.) 


-(2-7)/(7-1) 


K'/“») (1 + %2) aaci + It 


(22) 


Moreover, by applying the same iype of analysis as described in refer- 
ence 2 to the seventh-degree velocity profile, the following expression 
for A(§) (denoted by >tg) can be derived (details are given in the 
appendix) : 


\{l) =^c 

2Fis 






(2/jis)«Pls [(m)/ (7-1 ) - (91s/Fis)1 

(^iN pi/T„)'- ^ 


(23) 


8lt should be noted that such a seventh-degree velocity profile is 
used here only for the purpose of determining the separation point and 
that otherwise the sixth-degree velocity profiles of reference 2 (see 
also eq. (Ij-)) should be used. 

9ihe use of this approximate value greatly simplifies the calcula- 
tions without appreciably afTecting the accuracy. This is due to the 
fact that equation (lO) should ordinarily be a fairly good approx im ation 
for bj (cf., e.g., table Hi), while, except in cases of extreme 
cooling (Gi close to zero) , the b3_ term in equation (20) will be 
relatively small. 
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where Fls and qp^ are constants given hy 

Fjs = 0.1159 + 0.002525a2s - 0.001l4-54a2g^ - 0. 0000572 (hiagg/G^)^ - 
0.0005T^('bia2s2/Gi) + 0.000887(hia2g/Gi) 

► 

cf^ = O. 25 G 1 + O.Ol^57 + 0.0758b2^ + 0.0548a2s - 0.00291a2g^ + 

0.00775 (bia2s/Gi) - 0.001ll)-7('bia28V'^) ~ 0*00011^5 (biagg/Gi)^ 

( 21 ^) 


ant^ h^ and a2g are 8ln5)le functions of according to equa- 
tions (10) and (20) . The constants ^Is and are functions of 

Gi and are shown in figure if. 

For any given reference Mach number and uniform wall-teinperature 

ratio Gi, the separation point in a region of given adverse pressure 
gradient, as specified hy uq^/UoqC &) f will he the station | at which the 
ri^t sides of equations ( 22 ) and (23) are equal. Baus, it is necessary, 
in general, only to plot A versus in the anticipated vicinity of 
separation, in accordance with hoth equations ( 22 ) and (25) } and to 
determine the point of intersection of these two cxirves. !Ehe separation 
point will evidently he independent of C, so that for the purpose of 
determining the separation point one my set C = 1 . In ilie case of a 
region of favorable pressure gradient steirting at the leading edge 
followed hy a region of unfavorahle pressiire gradient, equation (25) 
should he modified to equation (A 9 ) of the appendix, "based on the 
assumption that A is continuous at the point of discontinuity of the 
velocity gradient Uj^ ’ . 


General Tmp li cations 

With respect to the effect of wall temperature on the separation 
point, it should he noted, first of all, that, as remarked in the pre- 
ceding paragraph, the value of the teinperature-viscoslty factor C has 
no Influence on the separation point. Thus, the temperatiJre-vlscosity 
relation, as it is incorporated in C, does not have the significant 
effect on the separation point that it has on the skin-friction eind 
heat-transfer coefficients- The effect of wall tenqjerature on the 
separation point ^d .11 thxis arise only frcan its dy n a mi cal, or pressure- 
gradient, effect (cf . the section "Skin-Friction and Heat-Transfer 
Characteristics") . 
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For a fixed velocity distribution U2yUj^(0 outside of the houndary 
layer and a fixed Mach nianber the equations developed here will 

imply that diminishing the wall temperature, that is, diminishing Gi, 
will have a favorable effect on separation by moving the separation point 
downstream. This can be seen particularly from equation (22), according 
to which the required value of X for separation (Xggp) will increase 
as Gi is diminished. This delay of separation caused by cooling of 
the wall is an illustration of the general tendency, discussed previously, 
of a decrease in wall temperature to diminish the direct Influence of a 
pressure gradient (in the present case, an adverse pressure gradient) . 

The effect of Mach nttmber on the separation point for a fixed value 
of U]yUoo(5) and either a fixed value of or a fixed value of 

To/^lio csinnot be predicted quite so readily from the equations as the 
foregoing effect of wall temperature. Mach number effects are contained 
in the Ti/T„ terms in eqioations (22) and (25) , as well as in the 
term of equation (22) . These terms tend to cause a decrease in both 
Xggp and Xg wiili increase in Mach number. For a fixed value of G^, 

that is, a fixed ratio of actual wall temperature to equilibriTmi wall 
tenqjerature for an insulated wall, nvnnerical exanq)les for the case 
ui/Uoo = 1 - S have indicated that the effect of Mach number on the 
required value of X, that is, Xggp (rather t h a n on Xq), is the pre- 
dominant effect, so that under such conditions an increase in Mach 
mmiber moves the separation point forward and thias enhances separation. 

A well-known case in this respect is that of zero heat transfer, where 
Gi = 1 (e.g., refs. 4 and I8), However, according to additional resvilts 

for the case uqyu,^ = 1 - |, lowering the fixed value of G^ tends to 
diminish somewhat this adverse effect of Mach nimiber (ref. l) . 

If, Instead of considering the ratio Gq as fixed, the ratio 
of wall temperature to the reference, or free-stream, tenperature is held 
fixed, then the effect of Mach number on the separation point for a given 
veloclly distribution ui/iaj„(5) is altered. This is due to the fact 
that, if To/^4o is fixed, then, as seen from eqmtion ( I 5 ) , G^ varies 
with Mach number, and, in particular, G]_ decreases as increases. 

This effect has a tendency, according to equation (22), to Increase the 
value of X required for separation as is Increased. This, in 

turn, tends to move the separation point downstream, and the numerical 
example to be subsequently given here Indicates that the net effect, at 
least in the case = 1 - |, of increasing -the Mach number as 

the wall-tenperature ratio Tq/T,, is fixed is to delay separation. This 
is hence in contrast with the results obtained for fixed values of G^^ 
that is, fixed values of Tq/Tq. 
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Numerical Exanqjle 

To illustrate quantitatively the application of the equations 
developed here, as well as the foregoing conclusions, a numerical 
example based on the ease 


^l/Uoo = 1-1 (25) 

will now he given in detail. Equation (25) represents the simplest type 
of an adverse pressure gradient. ^ 


First, the effect of wall tenperature on the separation point for 
a fixed Ifech number will be calculated. For this purpose, it will be 
assumed for simplicity that ly^j = 0 (so that T^yo^o = 1) while the 
parameter Gq (Gq = Tq/Too if = O) varies. In this case, equating 
the right sides of equations (22) and (25) leads to the following sinple 
expression for the location of the separation point as a function of the 
wall tenperature (I = Isep) • 


I 


sep 


= 1 - 


^ ^ ^?ls \ (“^Is/^'Pls) 

. + V 


(26) 


where arifl are given in terms of G^^ by equations (24), (lO), 

and (20) or by figure 4. A result qualitatively similar to equation (26) 
was obtained by the use of fourth-degree profiles in reference 1, but 
becaxise of the present use of a seventh-degree velocity profile the 
resvilts based on equation (26) are considerably more accurate quantita- 
tively. The separation point according to equation (26) is plotted in 
figirre 5 for a range of valiies of 3?o/^oo from 0.5 to 2.0, and the effect 

of cooling of the wall in delaying separation can be clearly seen here. 


In order to determine the effect of Mach nxmiber for a given ratio 
of wall tenperature to free-stream tenperature T^yi^, the sei)aratlon 
point has been calculated for a range of Mach numbers from 0 to 5-51 

with a fixed value of Tq/U^,,,, namely, = 2. The values of G^^ 

therefore, range from = 2 (for = 0, wall heated) to Gi = 0.3 
(for = 5.31, wan cooled). The calculations can be performed by 
observing that by virtue of equations (n) and (12) equation (22) with 
7 = 1.4 can be written as 


lOif ui/Uj^ = 1 - k|, where k is a positive constant, then, as 
shown in, for example, reference 4, the results based on equation (25) 
remain valid, with | replaced by k| and A replaced by kA. 
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Agep = -105c 


(T]yg^)-o-^ 

\k + U(To/Tj + 


(27) 


Equation (27) is in a form convenient for calculations involving a fixed 
value of Iq/Too- For each specified Mach number ^ep ®' 

function of | can thus he readiiy calculated, ■while AgCl) can be cal- 
culated from equations (23) and (2ll-) (or fig, 4 ) by numerical integration. 
The ■value of | at which Tvg = ^ep then the separation point. The 
results of such a calculation for -tiie case denoted hy equation (25) are 
shovm in figure 6. Fcr comparison, ■fche separation point ■versus Mach 
nuMjer for zero heat -transfer is also included in the figure. The 
results clearly indlca-te -fche favorable effect of Mach number (for -tlie 
fixed velocl-fcy dls-tribution of equation (25)) on ■tiie separation point 
for -fche fixed ■value of Tq/I^jo, in contrast with -fche unfavorable effect 
of Mach number at zero heat "transfer (when G]_ is fixed at uni-ty and 
hence Tq/Q^ increases ■wi-th Mach number: cf. eq. (15)). 


S^tagnation Flow Followed by Adverse Pressure Gradient 

It has already been indlca^ted ho\'T "the separa-tion point can, in 
general, be calcula-bed in cases of a favorable pressure gradient followed 
by an adverse pressure gradient. It may be of interest in -this connec- 
tion to investiga^te the following question: Under what conditions will 

■fche boundary layer in the favorable-gradient region develop to a suffi- 
cient extent so ■fchat laminar separation ■will occxir almost Immediately at 
■fche point where -fche unfavorable gradient s-fcarts? In particular, to ■what 
extent does -fche ■wall -fceiirpera^fcure affect such conditions? From a practi- 
cal point of ■view, -fchis question appears equi-valent to the question of 
when -fche flow ivlll separa-fce at -fche point of mini mum pressure outside of 
the bo\mdary layer. It is well known -that, at least in Incoanpresslble 
flow wi-fchout heat ■transfer, 1 am1 nar flow usually ■fcends to separate 
shortly downs'fcream of such a point if the pressure gradient is contin- 
uous. In case -fche favorable pressxire gradient can be represen-fced by a 
s-fcagnatlon flow, -fchat is, by a flow outside of -fche boundary layer of -fche 
form 


^l/H» = (28) 

■where K]^ is a positive cons-fcant, and ■fche Mach numbers in this region 
are sufficiently low so -fchat -fcheir effect in -fchis region can be neg- 
lected, it -will be seen -fchat -fche foregoing questions C£in be answered 
in a particularly sin^ile and fairly in-fceresting manner. 



20 


mCA OIN 5296 


let the subscript » now denote the point where' the adverse 
pressure gradient starts (here presumed abruptly, that is, discontinu- 
ously^) after the favorable pressure gradient, and let Kg be pro- 
portional to the mEignitude of the negative velocity gradient in the 
adverse-pressure-gradient region at this point, that is, 

K2=-(ui‘/u„) (29) 

S> =oo4- 

Then, according to equation ( 22 ) , separation will, in general, occur 
immeiiately at this point if X (denoted here as >^av) this point 
as calculated frcm the flow in the favorable pressvire gradient satisfies 
the relation^ 


>fav 2 —7 ^ (30) 

K2(i + tt.2) UCSi + It 

If, in particulEir, the region of favorable pressiire gradient is a stag- 
nation flow characterized by equation (28) , then, applying the results 
of reference 19 and identifying X with the quantity there. 


Vav = ^(Gl)/Ki (51) 

where f(Gi) is a function of which can be found either from 
fig\ire 1 of reference 16 or "by solving the algebraic equations (50) and 
(31) there. 15 Suibstitution for Xf^^ into relation (30) yields the 
following condition for immediate separation after the stagnation-flow 
region; 



^^■^lis is, of course, an idealization, since in actuality the 
pressure gradient will not be discontinuous. However, this idealiza- 
tion ml^t also be regarded as an approximation for a rapidly changing 
pressure gradient and serves to furnish at least a qualitative answer 
to the foregoing questions. 

^It is permitted to put C = 1 here since it has already been 
seen that C does not affect the location of the separation point. 

^ f(Gi) is essentially v.-^ for cp = 0 in reference 19. The 
quantity h there is equivalent to the quantity G^ here. 
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wliere 


F(Gi) 


105 

(11% + 4)f(Gi) 


(53) 


In figure 7 F(Gi) is plotted against G]^ for G^ = 0 to 1. As seen 
from figure J, F(Gi) varies from 2.40 to 1.00 in this range and 
increases as Gi d im inishes, that is, as the -wall is cooled. 

From relation (32) it is seen that the condition, in the present 
case, for separation at the immediate start of the adverse pressxire gra- 
dient depends in a sinple manner on the magnitude of the ratio of the 
adverse velocity gradient K2 to the favorahle velocity gradient %. 
Cooling of the wall is seen, once again, to have a tendency to prevent 
separation, since it increases the minimum required value of Kg/Ki. 

The Mach number at the beginning of the adverse-pressure-gradient region 
is, however, seen to have an unfavorable effect. 


STABILITY CHAEACTEIRISTICS 


It is of interest to investigate the effect of wall ten^erature, 

Mach number, and pressure gradient on the stability characteristics of 
the laminar boundary layer. For this purpose, two "types of calcula-blons 
will be made for "the srq)er sonic flow over "the thin airfoil (fig. l) on 
which the nimierlcal exsraple in "the section "Skin-Friction and Heat- 
Transfer Characteristics" was based. First, "the wall -tempera-ture 
required to s"tabillze "the laminar boundary layer conple"tely, "that is, 
for infinite mlniimnn critical Reynolds number, will be calcula-ted at 
"two given s"tatlons along "the flow for "various tfech numbers. Second, 

"the minimum critical Reynolds number at a gi"ven station will be calcu- 
la"ted as a function of "the wall tenperature. 

The me"thod of calculation is based on "the s"tabill"ty crl"teria de"vel- 
oped ty T.in and lees (refs. 5 and 6) wi"th certain modifications pre- 
sented in an unpublished paper entitled "Calculation of S"tabili"ty of 
Constant-Pressure Boundary Layers on Iso"thermal Surfaces Wi"th an In"tegral- 

Method Mean-Flow Solution" by Martin Blocm.^5 Explicit de"tadls on the 
method of calculation can also be foimd in reference 4. The stabili"ty 
criteria as developed by Lin and Lees are based on the amplification or 
decay of small disturbances, nnd the minimum critical Reynolds number 
thereby obtained is the minimum Reynolds mmiber required for "the 

^^Because of the neglectlon of Mach nunber effects in the s-tagnation- 
flow region, equations (32) and (33) may not be quantitatively "valid for 
hi^ "values of M^. 

15 

See also references 9 to 11. 
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possildllty that at least certain types (depending on the wave lengths) 
of disturhances will he amplified. Thus, the existence of a Reynolds 
number exceeding the minimum critical Reynolds number is a necessary 
condition for instahillty of the laminar boundary layer. However, 
transition to turbulence, which appears to depend on the magnitude of 
the aji^lified disturbances, will usually occur at a higher Reynolds 
nimiber and hence at a point downstream of the point where the actual 
Reynolds number Is equal to the ml ni Tmnn critical Reynolds number. Never- 
theless, it can probably be qualitatively concluded that the higher the 
minimum critical Reynolds number, the more stable the flow and the less 
the tendency for transition. (Further details on such questions can be 
found in references 5^ 6, and 20,) 

Figure 8 shows the wall-temperature ratio Tq/Tij and figure 9> "the 
ratio To/T]_ versuis the I&ch number required for ccmplete stabili- 

zation of the flow at two different stations along the airfoil of fig- 
ure 1. The resTilts for a flat plate (zero pressture gradient) are also 
included for comparison. Figure 9 clearly indicates, from one viewpoint, 
the stabilizing influence of the negative pressure gradient here, since 
the required maxiimjm values of Tq/T]_ for infinite nvinimum critical 
Reynolds number are greater for the airfoil than for the flat plate > 
hence, less cooling, relative to the local temperature T^^ outside of 
the boundary layer, would be required for the airfoil than for the plate. 
Figure 9 indicates that at the hl^er Mach niimbers the stabilizing 

effect of the favorable pressure gradient appears to be relatively 
diminished . I 6 

These results on the effect of the pressttre gradient are qualita- , 
tively in accordance with the conclusions of reference 12 (based on zero 
heat transfer) . For the special case of M^, = 1, it may be observed, 
in this connection, that, althou^ infinite cooling (Tq = 0) would be 
required to stabilize the boundary layer over a flat plate ccopletely, 
only a finite degree of cooling (Tq > 0) would suffice to stabilize the 
flow over the present airfoil. This is due to the fact that the local 
Mach number increases along the flow over the airfoil, so that when 
= 1 the local Mach nuniber Mp will exceed unity at the stations 
along the airfoil dOTnastream of the leading edge. 

It is significant to note that the c\irves for the required maximum 
values of To/Tb (fig- 8) for the airfoil and for the flat plate cross 
one another at certain Mach numbers This indicates that, for a 

fixed reference temperatm:e T^ outside of the boundary layer immedi- 
ately behind the shock wave at the leading edge, there are Mach 


i6lhis mi^t be due physically, at least in part, to the fact that 
for the higher tfeich numbers the magnitudes of the favorable-velocity- 

gradient ratio Ui'/Ub according to equations ( 17 ) and ( 18 ) as well as 
of the pressure-gradient parameter a^ (table IV) are decreased. 
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nimibers (for exar^jlLe, > 2-5 at |’ = 0.8) for -wMch a lower 

(uniform) wall tenqperature would le required on the airfoil than on the 
flat plate to stahlUze the flow ccaiipletely. This seemingly paradoxical 
result (in view of the favorable effect of the negative pressure gradient 
described In the preceding paragraph) Is due to the fact that the local 
temperature Tq outside of the boundary layer over the airfoil d^.mlni shes 
along the flow and Is therefore less than Tjj, partlculeirly for the higher 
Mach numbers (cf. eq. (U), with subscript » replaced by b) . Thus, 
from the viewpoint of ccraplete stabilization, the net effect of the nega- 
tive pressmre gradient at higher Mach numbers Is Tinfavorable with respect 
to the required temperature ratios To/Tb but favorable with respect to 
the required temperature ratios ^ sense of figures 8 and 9* 

Figure 10 and table VI show the minimum critical Reynolds nxmiber Rb,cr 

for the boundary layer at a given station along the airfoil and at two 
different tfeich numbers for both various wall -reference temperature 
ratios To/Tb and various wal.l -equilibrium temperature ratios G^. The 
stabilizing effect of cooling of the wall Is clearly Indicated here, since 
the TTi1nlTTmm critical Reynolds number Is seen to Increase as the wall- 
tenperatvure ratios are diminished. Moreover, by ccmparlson with the results 
for flow over a flat plate, the stabilizing effect of the negative pressure 
gradient (by Increasing the min im um critical Reynolds number) Is also seen 
In figure 10. The destabilizing Influence of a positive pressure gradient 
Is, for zero heat transfer. Illustrated by an example In reference k. 

It Is Interesting to observe the effect of Mach nvraber on the sta- 
bility of the .laminar boundary layer. From figure 10, It Is seen that 
for a fixed ratio of wall tenperatxnre to equlUbrlima temperature 

an Increase of Mach number from I .5 to 2.0 destabilizes the boundary layer 
both over a flat plate and over the airfoil. This effect, for the limited 
Mach nimber range considered. Is seen. In fact, to be enhanced by the 
negative pressure gradient here. If, now. Instead of a fixed value of 
G^, a fixed ratio of wall temperature to reference temperatinre To/Tb Is 
considered, the effect of Mach nxmber,' for the limited range under con- 
sideration, Is changed. For flow without a pressure gradient, an 
Increase of Mach number Is now seen from figure 10 to have a stabilizing 
effect, especially at the lower wall temperatures. For the flow over the 
airfoil, however, figure 10 (cf. also table Vl(a)) now Indicates that an 
Increase of Mach number has a stabilizing effect only at wall temperatures 
close to the critical wall temperature (l.e., for Infinite mlnimmjm criti- 
cal Reynolds n^mber) and that for (fixed) hl^er wall- temperature ratios 
of Tq/^ ^ Increase of Mach number has a clear destabilizing effect, 
similar to the case of fixed valxies of G^. Thizs, for the low Mach nimber 
range treated here, the negative pressure gradient over the airfoil con- 
sidered here modifies the effect of Mach number on laminar stability for 
a fixed ratio of wall temperature to reference free-stream temperature. 
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The criticELl EeynoMs numbers in figure 10 are "based on conditions 
immediately "behind the leading-edge shock wave on the airfoil. By 
comparing these Reynolds numbers mth those for a flat plate = l) , 

the effect of the actual sh.ock-T'ra.ve on the sta"bility of the "boundary 
layer is essentially eliminated, so that only the effect of the velocity 
(or, equivalently, pressvire) gradient eiddition to the effect 

of wan tengjerature) is included here- From a practical viewpoint, it 
may also "be of interest to detennine the effect on "boundary-layer sta- 
bility as the flight speed of the stqjersonic airfoil is increased. For 
this purpose, the minlrmim critical free-stream Reynolds number Roofer 
versus the free-stream Mach number has "been calculated, and the 

results are shown in figure 11 and table "71. The results are seen to be 
quite similar to those based on conditions immediately behind the leading- 
edge shock wave in the airfoil. 


CONCLUSIONS 


"Under the assunq)tions of a Prandtl number of 1 and a linear viscosity- 
temperature relation in conjimction with Sutherland’s equation, the 
follOTiTlng conclusions can be stated frean the present investigation of the 
ccanpresslble boundary layer in a pressure gradient over a airrface at a 
given uniform wall teng)erature. (interaction between the boundary layer 
and the external stream has not been considered.) 

1. The effect of wall teniperature on the skin-friction and heat- 
transfer coefficients arises from the pressure gradient and (independently 
of the press^Ire gradient), from the factor C in the viscoBlty-"ben 5 )erature 

relation. (C = (Tq/To,) + S)/(Tq + S) where Tq 1b the temperature 
at the wall, is the tenq>erature at a point outside the boundary layer, 

and S is the Sutherland constant.) In regard to the pres sure -gradient 
effect, cooling of the wall tends to diminish the Nusselt number and, 
especially, the skin friction in a favorable (negative) pressure gradient 
and to increase the coefficients in an adverse (positive) pressure gra- 
dient. In regard to the temperature-viscosity effect, lowering the wall- 
tenperature ratio Tq/Tm "will ordinarily tend to increase both the skin 
friction and Nusselt number. 

2. Cooling of the wall tends, in general, to dimlnlBh the direct 
effect of a pressure gradient. A particularly clear exanrple of this is 
the delay of separation in an adverse pressTire gradient by cooling of 
the wall. 

5- A simple and ordinarily siofficiently accurate method of deter- 
mining the separation point in a given subsonic or supersonic adverse 
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pressiire greidient over a -wall at any specified uniform wall ten 5 )erature 
has "been developed here. 

It. [Che results of a numerical example for a fixed linearly decreasing 
velocity outside of the houndeiry layer indicate, in addition to the 
delaying of separation hy cooling of the wall, that, for a fixed ratio Gi 
of wall temperature to equilihriim adiahatlc wall temperature, an increase 
of free-stream tfech number moves the separation point upstream; while, 
for a fixed ratio of wall tenperatiire to free-stream temperature, 

an increase of tfech number has, in general, a less unfavorable effect and 
in this case actually moves the separation point downstream. 

5- Numerical examples for the supersonic flow over a thin airfoil 
indicate in detail the stabilizing nature of a negative pressure gradient 
and of cooling of the wall on the laminar boundary layer. The examples 
also indicate that the pressure gradient here modifies the effect of Mach 
nimnber on laminar stability for a fixed ratio of wall temperature to 
reference free-stream temperature. 

Polytechnic Institute of Brooklyn, 

Brooklyn, N. Y., April 1-95^ • 
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APPENDIX 

DETSEMINAfEION OF SEPAEATION POINT 


Tftien the assungationa that the Prandtl number is 1 and that the 
coefficient of viscosity is proportional to the absolute temperature 
(eq. (1)) are made, differentiation of the momentum partial differential 
equation with respect to t yields the following relation at the wall 
(cf. eq. (25) of ref. 2 ), where u = v = 0 : 

/ 

Ut5jo 5?o\^t ^t2/o 



Moreover, the energy partial differential equation and differentiation 
of this equation with respect to t yield the following relations for 
a uniform wall temperature (cf . eqs. (26) and (27) of ref. 2 ) : 

= 0 

» 

(^H/St5)o = 0 

where 

H s (u 2/2) + CpT 

From equations (A 2 ) and (A 5 ) it follows that at the separation 
point, that is, where 

(Su/Bt)Q = 0 (a 4 ) 

the relation 

= 0 (A 5 ) 


(A 2 ) 


(A3) 


holds. Differentiating the momentum paxtial differential equation twice 
with respect to t and taking values at the Bepjaration point, it is 
found, with the use of equations (A 1 ) to (A 5 ) , that 

)o = 0 (a 6 ) 

which is in accord with equation (19) of the main text. 
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The seventh-degree velocity profile in t satisfying condition (a6) 
in addition to the boundary conditions (viz., eqs. (22) to (28) of ref. 2) 
satisfied in the present analysis hy the sixth-degree profiles is 

= (4^ - “ 1 '"’^) + + t 2 - |r5 + 

+ t5 - 5t5 + 

where a£ is given hy eqmtion (6), while 

83 = (82^l)/(5Gij (a8) 

The condition (8u/8t ) q = 0 will, according to equations (A7) anfl (a8) , 
lead to eqmtion (20) of the main text. 

By inserting the profile (A7) , in conjunction with equation (a8) , 
into the momentum integral-differential equation (U) of reference 2 
and assuming, as in reference 2, that the and hj^ terms in Fj 

and F2 (defined in ref. 2) may he replaced hy constant values, an 
ordinary differential equation of the same form as equation (9) of the 
main text is obtained, except that ^2. ^nd (now written as F]_g 

and q^) are now given hy equations (24) of the main text, while the 
factor 2 on the ri^t side of equation (9) is replaced hy 7/4. Com- 
parison, accordingly, with the solution of equation (9) (eq. (7)), leads 
to equation (25) of the main text. Since the chief purpose is here the 
location of the separation point, the constant value of (as in 

ref. 4) is now chosen as that at the separation point and hence as that 
given hy equation (20) of the main text. 

In case the region of adverse pressure gradient starts at some 
point 5 = 5a <3.ownstream of the leading edge, then equation (25) can 
still he applied directly in calcxilating the separation point. Greater 
accuracy, however, mi^t he obtained in such a case hy applying equa- 
tion (23) only for the region of adverse pressure gradient. For this 
purpose, equation (23) must he modified to satisfy the boundary condi- 
tion A = at i = Thus, 

2 27-1 ^is 

(I)|=g^ Aa + I C • (Ti/!B„) ^ d| 

\ (A9) 
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wtiere 

2 7+1 ^ 

I = (V“J ^ ( V’i.) • 

and where \ can he obtained as the value yf A at I = la "based on 
equation (?) for the region (O < | < la) of favorable pressure gradient. 
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TABLE I 

VALUES OF TEMPERAOJURE-VISCOSITY FACTOR C 
VERSUS 

MACH NUMBER AND WALL-TEMPERATURE RATIO 


% 

^ 1:5 1 

3.0. 

Gl 

0.3 

0.5 

0.7 

0.3 

0.5 

0.7 

C 

1.097 

1.057 

0.997 

1.033 

0.923 

0.834 


TABLE H 

VALUES OF yc ALONG AIRFOIL OF FIGURE 1 


1 ' 


VC at 

of - 



1.5 

5.0 

Gi = 0.3 

Gi = 0.5 

II 

0 

Gi = 0.5 

Gi = 0.5 

Gi= 0.7 

0 

0 

0 

0 

0 

0 

0 

.1 

5.73 

3.7J+ 

5.76 

5.92 

5 . 9 l^ 

5.95 

.2 

7.49 

7.49 

7.51 

8.27 

8.28 

8.28 

.3 

11.32 

11.28 

11.26 

15.07 

15.05 

13.01 

.4 

15.22 

15.11 

15.02 

18.54 

18.24 

18.14 

.5 

19.23 

19.01 

18.85 

24.12 

23.91 

25.70 

.6 

23.35 

25.01 

22.72 

30.41 

50.08 

29.71 

.7 

27.56 

27.09 

26.67 

37.29 

56.75 

36.18 

.8 

31.86 

31.22 

50.67 

44.50 

45.70 

45.00 

.9 

36.41 

35.56 

54.82 

52.50 

51.29 

50.32 

1.0 

41.01 

39.95 

59.01 

60.49 

59.18 

57.90 
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TABLE HI 

VALDES OF 131/(1 - Gi) ALCMG AIRFOIL OF FIGURE 1 


I* 



./(^ " ®l) 

at of - 


1.5 

5.0 

6 

II 

Gi = 0.5 

Gi = 0.7 

, 5 ^ 

II 

0 

Gi 

Gi = 0.5 

Gi = 0.7 

0 

2.000 

2.000 

2.000 

2.000 

2.000 

2.000 

.1 

1.978 

1-975 

1.-969 

1.980 

1.976 

1.974 

.2 

1-959 

1.955 

1-945 

1-963 

1.958 

1.949 

.5 

1.9Mj- 

1.935 

1.927 

1.947 

1.939 

1.929 

.k 

1.932 

1.921 

1.910 

1.935 

1.924 

1.913 

.5 

1.923 

1.910 

1.897 

1.924 

1.915 

1.900 

.6 

1.914 

1.901 

1.889 

1.916 

1.904 

1.900 

•7 

1.908 

1.895 

1.882 

1.910 

1.898 

1.883 

.8 

1.904 

1.890 

1.879 

1.905 

1.8^ 

1.879 

• 9 

1.901 

1.887 

1.877 

1.902 

1.891 

1-879 

1.0 

1.898 

1.886 

1.876 

1.900 

1.890 

1.876 


TABLE IV 

VALDES OF -E 2 ALONG AIRFOIL OF FIGDRE 1 





-a 2 at 

of - 



1' 

1.5 

5.0 


1^? 

II 

0 

Gi = 0.5 

Gi = 0.7 

Gi = 0.5 

Gi = 0.5 

Gi = 0.7 

0 

0 

0 

0 

0 

0 

0 

.1 

.1205 

.2017 

.2840 

.1077 

.1802 

.2536 

.2 

.2515 

.5859 

.5414 

.2128 

.5548 

.4972 

-3 

.5551 

-5567 

.7777 

.3341 

-5219 

-7296 

.4 

.4509 

.7130 

.9920 

.4115 

.6816 

.9492 

.5 

.5202 

.8574 

1.1890 

.5029 

.8506 

1.1527 

.6 

.6044 

.9928 

1.5721 

.5881 

.9696 

1.5408 

.7 

.6^98 

1.1140 

1.5555 

.6669 

1.0949 

1.5698 

.8 

.7513 

1.2271 

1.6877 

-7531 

1.2002 

1-6551 

-9 

.8175 

1.3305 

1.8259 

.7899 

1.2911 

1.7730 

1.0 

.8764 

1.4251 

1.9454 

.8591 

1.5565 

1.8579 
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TABLE ? 

SEPARATION POINT CALCUIAIED AS A FUNCTION OF 
MACH NUMBER FOR ZERO HEAT TRANSFER 
(Gi = 1) AND ux/Uoo = 1-^1 



^sep far of - 

0 

1 

5 

10 

Based on eqs. (22) to (2^) 

0.122 

0.115 

0.0768 

0.025 

Baaed on metBod of 
Stewaxtson (ref. I5) 

0.120 

o.uo 

0.077 

0.024 
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TABI£ VI 


MUJIMCM CRITICAL REYNOLDS NUMBERS OF LAMINAR 


BOUNDARY LAYER OVER AIRFOIL IN FIGURE 1 


AND OVER A FIAT PIATE 


(a) Over airfoil (|’ = 0 . 8 ) 


Ko 


1.662 


2 . 184 - 


To/Tb 


1.1 1.595 i-‘708 

1.0 1.4-50 1.553 

.93 1.34-8 l.Mi -3 

.91 1.320 1 . 4-13 

.90 1.305 1-397 

1.379 


1.0 1.8 1.94 

.9 1.62 1.74- 

.8 1 . 44 - 1.55 

.77 1.386 1.49 

.737 1-327 .1.43 


®b,cr 


2.69 X 10^ 
6.58 
48.6 
128.6 
261.3 


0.217 

.389 


00, cr 


2.59 X 10**- 
6.34 
46.86 
124.0 

251.9 


0.201 

.362 

1.653 

4.137 


(L) Over flat plate 




1.2 

1.0 

.95 

.90 

.75 

.72 

.7166 


1.00 

.90 

• 75 
.70 
.683 



^,cr 

cr 

1.740 

0.0804 X lo 4 

0.0804 X io 4 

1.450 

.282 

.282 

1.578 

.426 

.426 

1.305 

.641 

.641 

1.088 

16.68 

16.68 

1.044 

322.2 

322.2 

1.039 

00 

00 

1.8 

0.0506 

0.0506 

1.62 

.0905 

.0903 

1-55 

.481 

.481 

1.26 

1.648 

1.648 

1.229 

00 

00 
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Figure 2 .- Distribution of local Hfusselt number along thin supersonic 

airfoil and along flat plate. 
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Figure 8.- Wall -reference -ten^jerature ratios Tq/^ required for infinite 

mini imim critical Reynolds mnriber. 
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Figure 9.- Ratios of wall to local main-stream wall ten^jerature 
required for infinite minimum critical Reynolds number. 
































